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Abstract. We define a doubly stochastic operator on a finite dimen- 
sional simplex and study the limit behavior of the trajectories under 
doubly stochastic operators. We prove that except for certain points, 
the trajectory of a point, under the doubly stochastic operator, tends to 
the center of the simplex. 
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1. Introduction 
A quadratic stochastic operator is defined as: 

m 

where 

m 

x e 5 m_1 = {x = (xi,x 2 , ■■■ , x m ) € R m : Xi>0, Mi = l,m, V] X{ = 1} 

i=l 

and the coefficients Pij k satisfy the following conditions 

rn 

Pij,k = Pji,k > o, y~]pij,k = i- 
*!=i 

This kind of operators arise in many models of population genetics [5] 
and the main problem is to study the dynamics, that is the limit behavior 
of the trajectories of quadratic stochastic operators which is rather difficult. 
Note that a lot of papers are devoted to studying this kind of operators 
[2, 8, 10, 11]. In [2] the class of Volterra quadratic stochastic operators 
was defined and the limit behavior of the trajectories was studied using the 
Lyapunov's function and the theory of tournaments. Moreover, the class of 
dissipative quadratic stochastic operators was outlined in [8] and the limit 
behavior of the trajectories was fully studied in [10]. In [3], the class of 
doubly stochastic quadratic operators was outlined, but the study of limit 
behavior of the trajectories is wide open. However, some results have been 
obtained for doubly stochastic quadratic operators. Namely, the necessary 
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and sufficient conditions for doubly stochasticity of a quadratic operator and 
extreme points of the set of doubly stochastic quadratic operators were stud- 
ied in [4, 7]. Moreover, the ergodic theorem for doubly stochastic operators 
was given in [9]. Even though, the ergodic theorem is proved for doubly sto- 
chastic operators, it does not give any information about the limit behavior 
of the trajectories of doubly stochastic operators. Therefore, the study of 
the limit behavior of the trajectories of doubly stochastic operators is still 
open. 

The aim of this paper is to define an arbitrary (not necessarily quadratic) 
doubly stochastic operator and study the limit behavior of the trajectories 
of doubly stochastic operators. We give some necessary preliminaries in the 
next section. In section 3 we give the main result of the paper. 

2. Preliminaries 

In this section we give some definition from majorization theory and give 
the definition of doubly stochastic operator. 

We define the (m — 1)— dimensional simplex as follows. 

m 

S m - 1 = {x = (xi, x 2 , • • • , x m ) e R m : Xi > 0, Vi = T~m, J2 Xi = 1 ^ 

i=l 

The set 

intS m ' 1 = {x G S™' 1 :xi>0, i = T~^} 
is called the interior of the simplex. 

For any x = (x\ , X2 ■ ■ ■ x m ) G s m ~ l due to [6] , we define xj. = (x^ , X[ 2 ] , ■ ■ ■ X[. 
where x^j > X[ 2 ] > • • • > ar[ m ]- nonincreasing rearrangement of x. The point 
x^ is called the rearrangement of x by nonincreasing. Recall that for two 
elements x, y taken from the simplex S m ~ l we say that an element x ma- 
jorized by y (or y majorates x), and write x ~< y(or y >~ x) if the following 
holds 

k k 

i=i i=i 

for any k = l,m — 1. In fact, this definition is referred as a weak majoriza- 

m m 

tion [1], [6], the definition of majorization requires x [i\ = Yl VU] - However, 

i=i i=i 

since we consider points only from the simplex, we may drop this condition. 
We refer to monographs [1], [6] for more information concerning majoriza- 
tion. 

A matrix P = {Vij)ij=~rm ls called doubly stochastic (sometimes bistochas- 
tic), if 

Pij > 0, Vi,j = l,m 

m m 

= 1, Vj = l,m ^Pij = 1, Vi = l,m. 
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If the entries of a doubly stochastic matrix consist of only O's and l's, 
then the matrix is a permutation matrix. 

A linear map T : S m ~ l — > S m ~ l is said to be T— transform, if T = 
XI + (1 — A)P, where / is an identity matrix, P is a permutation matrix 
which is obtained by swapping only two rows of I, and < A < 1. 

Lemma 2.1. [6] For the concept of majorization and x,y G S m ~ 1 , the 

following assertions are equivalent. 

k k 

1) x <y, that is £ x$ < ]T y [{ \, k = l,m - 1. 

i=l i=l 

2) x = Py for some doubly stochastic matrix P. 

3) The vector x belongs to the convex hull of all m\ permutation vectors 
ofy- 

4) The vector x can be obtained by a finite compositions ofT— transforms 
of the vector y, that is, there exist T— transforms Ti,T2,- ■ ■ , such that 
x = TiT 2 • • ■ T k y. 

So from the above lemma it follows that doubly stochasticity of a matrix 
P is equivalent to Px -< x for all x G S m ~ l . Motivated by this, in [3], the 
definition of doubly stochastic operator was given as follows. 

A continuous stochastic operator V : S" 71-1 — > S" 71-1 is called doubly 
stochastic, if 

Vx -< x 

for all x G S* m_1 . Identity operator, permutation operators (that is the 
linear operators with permutation matrix) and T— transforms are all doubly 
stochastic. The operator V : S 2 —> S 2 given by 

Vx = x 2 + 2yz, 
Vy = y 2 + 2xz, 
Vz = z 2 + 2xy. 

is an example for doubly stochastic quadratic operator. Doubly stochasticity 
can be shown using the definition or can be obtained as the consequence of 
Theorem 2.6 of [4]. 

Let V be doubly stochastic operator and x° G S m ~ 1 . The sequence 
V(x°), V 2 (x°), ...} is called the trajectory starting atx°. Here V k+1 (x°) = 
V{V k {x )), k = 1, 2, ■ ■ ■ and we also put V°(x°) = x°. We denote by u(x°) 
the set of limit points of the trajectory and it is said to be the u— limit set 
of the trajectory. We say that an operator V satisfies the ergodic theorem 
, if the limit lim X+Vx ^ vyn — £. exists for any x G S" 71-1 . For a doubly 

n— too n 

stochastic operator on a finite dimensional simplex the oj— limits set of any 
initial point is finite, hence V satisfies the ergodic theorem [9]. However, 
the convergence of Cesaro mean of the trajectory of dissipative q.s.o. does 
not give any clue about the convergence of the trajectories. Therefore, this 
problem will be studied separately in the following section. Note that the 
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dissipative operator was defined as Vx >- x in [8]. The limit behavior of the 
trajectories of dissipative quadratic stochastic operators was studied in [10]. 

3. The limit behavior of trajectories 

In this section study the limit behavior of the trajectories of doubly sto- 
chastic operators. 

A continuous functional (p is said to be a Lyapunov function for the opera- 
tor V, if the limit lim ip(V n (x )) exists along the trajectory {x°, V(x°), V 2 (x°) 

The Lyapunov function is considered to be very useful in the study of limit 
behavior of discrete dynamical systems (see [2] for example) . 

Theorem 3.1. A continuous functional given by 

m 
i=l 

is a Lyapunov function for doubly stochastic operator V. Moreover, if x £ 
intS m ~ l and Vx is not the permutation of x, then ip(Vx) < <p{x). 

Proof. Use the fact (see [6], prop. F.l, page 88.) that the function ip(x) = 
Y^, XiXj satisfies the condition that if x -< y then i/j(x) > if)(y). Moreover, 

i<j 

if x,y £ intS m ~ l , and x is not the permutation of y, then ip(x) > 4>{y)- 
Taking into consideration ip(x) = 1 — 2ip(x), we obtain that x ~< y implies 
ip(x) < ip(y). Since V is doubly stochastic, then Vx -< x hence (p(Vx) -< 
(fi(x). Further, since the simplex is compact, then the continuous functional 
ip is bounded. Note that the sequence a n = ip(y n (x )), where V is doubly 
stochastic and x° <G S™ -1 , is monotone. Hence the limit lim (p(V n (x )) 

n— s-oo 

exists. 

Now we show the second part. If x G intS m ~ l , then according to Lemma 
2.1 it follows that Vx <G intS m ~ 1 . Since all components of Vx and x are 
positive and Vx is not the permutation of x, then by applying the above 
mentioned fact, we get ip(Vx) < tp(x) 

□ 

By the same way one can prove that any continuous convex (or concave) 
symmetric ( (p is symmetric, if ip(Px) = (p(x) for any permutation matrix 
P ) functional is a Lyapunov function. 

We define the following set 

Ay = {i£ S m ~ 1 : 3i € N, and permutation matrix P such that V l (x) = PV l 

Theorem 3.2. Let V be doubly stochastic operator and x° G intS' m ~ l \ Ay- 
Then the trajectory of x° under V tends to the center of the simplex, that 
is, to the point 
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Proof. We define the following functional: 

m 

According to the Theorem 3.1 ip(x) is a Lyapunov function. Let us put 

lim <p{V n {x )) = C. 

n— >oo 

Since <p(x) is a Lyapunov function and monotonically decreasing along any 
trajectory, then the set 

<f~ l (C) = {xe intS™' 1 \ A v : <p(x) = C} 

is invariant for the operator V. 

Note that the minimum value of the function ip(x) is at the point 

If C + h then for x G ™tS m - 1 \ Ay according to the 
second part of Theorem 3.1 we have (p(Vx) < ip(x). This contradicts to the 
fact that (/9 _1 (C) is invariant under the operator V. Therefore C = — and 
<£ _1 (C) = {(^, ^, ■ ■ ■ , ^)}, which means that the trajectory tends to the 
center of the simplex. 

□ 

Corollary 3.3. The theorem holds true for those points x° G S™ 1 " 1 for 
which there is a natural number r such that V r (x) € intS m ~ l \ Ay. 

In fact, the sequence in Theorem 3.2 need not necessarily to be the orbit 
of doubly stochastic operator. One can generalize the above theorem as 
follows. 

Corollary 3.4. Let {x^} n >i be the sequence from intS m 1 such that 
y x (n+i) _j x (n) an( ^ x (n+i) i s nQ ^ ^ e p ermu tation of x^ for any n > 1. 
Then x^ ->• ( — , — , • • • , — ) as n ->■ oo. 

Proof. The proof goes along with the same ways as in Theorem 3.2. □ 

Now we are going to study the set Ay. If V is a permutation operator, 
then Ay = S' m ~ l . So the above theorem does not make a sense. However, 
the study limit behavior of trajectory of permutation operators is quite easy 
since the permutation operators act as a cyclic shift. Surely, except for 
permutation operators, the set Ay is not the whole simplex. 

The point x° is called p— periodic, if there is a number p such that 
V p (x°) = x°, and V t (x°) / x° Vz = l,p — 1. If p = 1, we say that it is 
fixed. The point x° is called eventually p— periodic if there is a number r 
such that V r (x°) is p— periodic. We just say periodic (or eventually peri- 
odic) if the period is irrelevant. Obviously, any periodic point is eventually 
periodic. 

Lemma 3.5. Eventually periodic points of the operator V belong to Ay. 
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Proof. Let x be eventually periodic point for V. Then there exist r such 
that V rp (x°) = V r (x°) for some number p. Since V is doubly stochastic the 
relation 



V r (x°) = V rp (x°) ~< V r{ P- l \x Q ) <■■■< V 2r (x°) ~< 
■< V 2r ~ 1 (x°) V r+1 (x°) ■< V r (x°) 

implies that there is a permutation matrix such that V r+1 (x°) = PV r (x°), 
hence x° G Ay. 

□ 

Theorem 3.6. For doubly stochastic operator and any initial point x° G 
intS m ~ l , which is not eventually periodic, the trajectory {x°, V[x ), V 2 (x°), ...} 
tends to the center of the simplex. 

m 

Proof. The function tp(x) = x i * s monotonically decreasing along the 

i=i 

trajectory {x° ,V(x°),V 2 (x°), ...}. Therefore lim tp(V n (x )) = C and hence 

to(x°) G (y9 _1 (C) Vx° G S m-1 , which means that the trajectory of any x° G 
S" 71-1 is convergent. The trajectory {x°, V(x°), V 2 (x°), ...} of the point x° 
can be rewritten as 

r^,o p o o p o o ... p o o l p i l p i l ... p i l 2 . .1 

where i^, j G N, i = 0, 1, ■ ■ ■ are permutation matrices, fcj > 0, j G N, 
and each x' l+l is not the permutation of x l (otherwise we can add one more 
permutation matrix to the sequence P[, P^ ■ ■ ■ P£.) . Note that the maximum 
value of kj, (j G N) is m\ (because there are m\ permutation matrices of 
order m). If there is a number j such that PjX 1 = x l ,i = 0, 1, • • • then 
the point x° is eventually periodic. Since we exclude eventually periodic 
points, then kj < ml for any j G N. This means that there exists an infinite 
subsequence {x°, x 1 , x 2 , ■ ■ ■ } of the trajectory {x° ,V(x°),V 2 (x°), ...} such 
that x 1 ~< x i_1 , % G N and each x l+l is not the permutation of x\ Applying 
Corollary 3.4 we get 



x ->•(—,—,■■■,—) 
mm m 



as i — > oo. Consequently, 



fV) ->(-,-,-,-) 
mm m 



as z — ? co. 



□ 



Corollary 3.7. Under the doubly stochastic operator, the trajectory of any 
initial point from the interior of the simplex either converges to the center 
of the simplex, or to a finite cycle. 



THE DYNAMICS OF DOUBLY STOCHASTIC OPERATORS 



7 



Let us now consider some examples. In these particular examples, the 
initial point need not to be from the interior of the simplex. 
Example 1. Let V : S 2 -> S 2 be given by 

Vx = x 2 + 2yz, 
Vy = y 2 + 2xz, 
Vz = z 2 + 2xy. 

The operator V is doubly stochastic [4] and has 4 fixed points 

(i,o,o) 1 (o,i,o),(o,o,i),(H,i). 

The first three are repelling and the last is attracting. If we define the 
function ip : S 2 — >■ R as tp(x) = tp(x, y, z) = x 2 + y 2 + z 2 , then one can see 
that 

2ip(Vx) = 3ip(x) 2 - 2<p(x) + 1 = (3<p(x) - l)(<p(x) - 1) + 2tp(x). 
Moreover, 

{X + \ +Z? <x 2 + y 2 + z 2 <{x + y + z) 2 = l 

implies that the maximum value of <f(x) is 1 at fixed points (1, 0, 0), (0, 1, 0), (0, 0, 1) 
and minimum is | at the point (|, |, |). If (p(x°) ^ 1, | then 

2<p(V(x )) = {3<p{x°) - l)(<p{x°) - 1) + 2<p{x°) < 2<p{x°). 

This means that the trajectory of any initial point, except for fixed points 
(1, 0, 0), (0, 1, 0), (0, 0, 1), converges to the center of the simplex. 
Example 2. Let V : S 2 -»• S 2 be given by 

Vx = z 2 + xy + xz, 
Vy = y 2 + xy + yz, 
Vz = x 2 + xz + yz. 

This operator is doubly stochastic [4]. The set of fixed points is {(x, y, z) G 
S 2 : x = z}, so it has infinitely many fixed points. One can also see that 
all the elements of the set {(x,y, z) G S 2 : y = 0} are 2— periodic. Simple 
calculations show that the inequality 

(Vx) 2 + (Vy) 2 + (Vz) 2 < + y 2 + z 2 

can be reduced to 

-2y(x - z) 2 (x + z) <0. 

The equality holds when y = which implies that the point is periodic, or 
when x = z, that is when the point is fixed. In all other cases, the inequality 
is strict. This brings to the conclusion that except for periodic points the 
trajectory tends to the center of the simplex. 

The example 2 shows that it is possible to find periodic points of a doubly 
stochastic operator on the simplex with period greater than 1. Moreover, 
it is possible to find 1— periodic (that is fixed) points in the interior of the 
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simplex. More generally, the center of the simplex is always fixed for any 
doubly stochastic operator (see [9], Th.l). Surely, all points of the simplex 
are periodic for permutation operator. However, from Theorem 3.5 there 
naturally arises the question whether doubly stochastic operators (which are 
not permutation operators) possess periodic points with the period greater 
than one in the interior of the simplex. Unfortunately, we are not able to 
answer to this question. 
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